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1. Introduction 

Recently, several new Chern-Simons matter theories with a large amount of supersymme- 
try have been found. In Gaiotto and Witten found a class of theories with A/" = 4 
super symmetry. Closely related are the the theories of ABJM |^], with J\f = 6 SUSY, and 
BLG with M = 8. All of these theories are superconformal, and arise as the low energy 
effective actions on certain brane configurations in string theory and M-theory. 

Some of these theories are conjectured to be holographically dual to M-theory on AdSj 
orbifold backgrounds, in similarity to = 4 super Yang-Mills theory in four dimensions, 
which was conjectured to be dual to Type IIB string theory on AdS^ x S^. Recall that, in 
the latter case, certain supersymmetric Wilson loops are dual to fundamental strings. Thus, 
as a check of this duality, one can compute the expectation values of these Wilson loops 
and compare to calculations made in string theory. This is difficult to do perturbatively, 
as the perturbative region of one theory is the strongly coupled region of its dual. 
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However, supersymmetric operators are often easier to deal with than their less sym- 
metric counterparts, and in it was shown that this is indeed true of supersymmetric 
Wilson loops in A/" = 4 SYM theory. It was demonstrated, using localization, that finding 
the expectation value of these operators reduces to a calculation in a matrix model. This 
allows one to compute it much more efficiently at any coupling, and the result provides a 
non-trivial test of the duality. 

In this paper we seek to apply the methods of to the J\f = 2 supersymmetric Chern- 
Simons matter theories discussed above. Our main result is that the partition function 
for a supersymmetric Chern-Simons theory with gauge group G and chiral multiplets in 
representations Ri, Ri* , R2, R2* , ■■■ localizes to the following matrix integral: 

Z = J (ia e"^*'^^'^'^^detyid(sinh(7ra))det_Ri(sech(7ra))deti?2(sech(7ra))... (1.1) 
where the integration is over the Cartan of the Lie algebra of G, and: 

detRf{a)=l[fip{a)) (1.2) 

p 

where the product runs over the weights of the representation R (or in the case of the 
adjoint representation (i? = Ad), over the roots of the Lie algebra). The same technique 
applies to more general representations, although the resulting matrix models are more 
complicated. 

We also consider the following supersymmetric Wilson loop: 



W = ^jT^trR (^Vexp (It (M^x^ + (j\x\)^^ 



(1.3) 



where a is an auxiliary scalar in the vector multiplet. We find its expectation value is given 
by: 

<W >= ^^^^j^ j do, e"'^*''''^'^"''trR(e^™)detAd(sinh(7ra))detiji(sech(7ra))detH2(sech(7ra)). 

(1.4) 

One application of this result is to a trivial example of a supersymmetric Chern- 
Simons matter theory: Chern-Simons theory without matter, which can be written in a 
supersymmetric form Here we are able to use localization to recover some well known 
results on the reduction of the Chern-Simons partition function to a matrix model 0, 
as well as reproducing some very simple knot invariants Q . 

Another example is ABJM theory. This is conjectured to be dual to a certain orbifold 
background in M-theory, so it would be interesting to make non-perturbative calculations 
in this theory. Here we were able to reduce the path integral to a matrix model, although we 
were not able to compute the resulting matrix integrals exactly. However, evaluating them 
as a perturbative expansion in the 't Hooft coupling, we find agreement with a perturbative 
calculation done in field theory |^, 1^, 11 1, which provides a check of our result. It is possible 



that the matrix model could be solved exactly in the large N limit using a saddle point 
approximation |12], as we will briefly mention at the end of the paper. 
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2. Setup 

The class of theories we will be considering, N = 2 supersymmetric Chern-Simons theory 
with matter, are described in for Minkowski space. We will work in Euclidean space. 
In this section we briefly review these theories. 

We start with the M = 2 gauge multiplet. This consists of a gauge field A^, two real 
auxiliary scalars a and D, and an auxiliary fermion A, which is a 2-component complex 
spinor. This is just the dimensional reduction of the N = 1 vector multiplet in 4 dimensions, 
with a being the reduction of the fourth component of All fields are valued in the Lie 
algebra g of the gauge group G. 

The kinetic term we will use for the gauge multiplet is a supersymmetric Chern-Simons 
term. In flat Euclidean space, this is: 

S = j d^xTi (^^"'P (^A^d.Ap + '^ApA.A^j - A+A + 2Da^ (2.1) 

Here "Tr" denotes some invariant inner product on g. For example, for G = SU{N), we 
will take "Tr" to mean k/iir times the trace in the fundamental representation, where k is 
constrained by gauge invariance to be an integer. 

This action is invariant under the usual (euclideanized) Af = 2 vector multiplet trans- 
formations: 

SAf, = I (r/t^^A - At7^e) 
5a = -i (??tA + Ate) 

5D = § ivh^D^A) - {D,A^h>'e) - \ {ri^\\a\ - [At,a]e) (2.2) 
b\ = {-ll^^'F^u -D + i-ft'D^a) e 

at = r/t {\-f^'^F^, -D- ij^'Dpa) 

Here e and r] are 2-component complex spinors, in the fundamental representation of the 
spin group SU{2). We will take 7^ to be the Pauli matrices, which are hermitian, with 
Ifiu = ^b^l,7u] = ie-tivpY- Also, Dp_ = + i[A^, .] is the gauge covariant derivative. 

Note that, in contrast to the Minkowski space algebra where we would have e = r/, in 
Euclidean space e and t] are independent. Taking e = rj would reduce us to the M = 1 
algebra. This is because there is no reality condition on spinors in 3 dimensional Euclidean 
space, so the least amount of supersymmetry one can have is a single complex spinor. 
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To carry out the locahzation, we will work on a compact manifold rather than in flat 
space, as this makes the partition function well-defined. As the above action is conformal, 
we can transfer it to the unit 3-sphere, S^, without changing any of the quantities we are 
interested in computing. The Lagrangian simply acquires an overall measure factor of ^^g. 
In addition, we must modify the following supersymmetry transformations: 

6D ^ ... + i(V^7?t^MA _ At7'^V^e) 

5X ^ ... + f a7'^V^e (2.3) 

<5At ^ ... - f CTV^7?t7/^ 

where now all derivatives are covariant with respect to both the gauge field and the usual 
metric on S^. One can easily check that this leaves the action invariant for arbitrary spinors 
e and r]. When we add matter, however, it will turn out to be necessary that we take e and 
ry to be a Killing spinors, which means that they satisfy the following equation: 

V;.e = ^^e' (2.4) 

Here e' is an arbitrary spinor. Note that in d dimensions, this is actually d equations, 
one of which determines e', with the rest imposing conditions on e. We will give the explicit 
solutions to this equation on below. With e a Killing spinor, the above supersymmetries 
give a representation of the superconformal algebra, anticommuting with each other to 
conformal transformations. 

2.1 The Wilson Loop 

The operator we will be localizing is the following supersymmetric Wilson loop: 

W = ^jT^tr/j (vexp (^j^ dT {iA^x" + a\x\)^^ (2.5) 



This operator has been considered in Here x^{t) is the closed world- line of the 
Wilson loop, and '"P" denotes the usual path-ordering operator. The variation of this 
operator under the the supersymmetry (|2.2[) is: 



5W oc -^r/t {^^x^ + \±\) A + ^At (7^0;^ - |x|) e (2.6) 
For this to vanish for all A we must have the following two conditions: 

T]^ {if^X^" + \x\) =0 

(2.7) 

(7^x^-|i|)e =0 

Note we cannot take e = rj here, which is why we need to consider theories with at least 
M = 2 supersymmetry. 

Now we need to impose the condition that e and r] are Killing spinors. This will force 
us to consider only certain loops, which will turn out to be great circles on S"^. To see this, 
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we will need to determine the Killing spinors on S^. We start by picking a vielbein. It 
will be convenient to use the fact that is, as a manifold, the same as SU{2), so we can 
take a local orthonormal basis of left-invariant vector fields ef. In terms of these, the spin 
connection is simply: 

= ^ijk^^ (2.8) 

Thus the spinor covariant derivative is: 

(2.9) 

We can immediately see a few solutions to the Killing spinor equation. Namely, take 
the components of e in this basis to be constant, in which case: 



V;.6 = |e^7fce = (2.10) 



This gives two of the Killing spinors. There are two more solutions, which can be seen 
most easily using a right invariant vielbein, and which satisfy: 

V^e = -'-j^e (2.11) 

Note that in these cases, e' is proportional to e. This is not true of a general Killing 
spinor (eg, take a linear combination of the above spinors), although in spaces of constant 
curvature it is always possible to form a basis of the space of Killing spinors with such 



special ones |14, 15 



Now let us impose the condition that e preserves the Wilson loop. If we pick r to be 
the arc length, we find that e must satisfy: 



i^^x^" - l)e = (2.12) 

This is only possible if 7^x^ is constant, which means must be some fixed linear 
combination of the e*, so we may as well pick our loop so that is parallel to one of them, 
say e^. The integral curves of these vector fields are great circles, so the Wilson loop must 
be a great circle to preserve any supersymmetry. Then this equation becomes: 

(73 - l)e = (2.13) 

So this restricts us to only one of the two left-invariant Killing spinors (there is also a 
right handed one that preserves it). We could also pick r] to be one of these Killing spinors. 
Thus this Wilson loop preserves half of the super symmetries. 

Conversely, given a Killing spinor e, and taking r] = 0, there is a family of great 
circles such that the Wilson loops along these circles are preserved by the corresponding 
supersymmetry. These are just the integral curves of the vector field e^7^e, which, since 
this vector field is left-invariant, form a Hopf fibration. As a result, one could use the 
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localization described here to compute the expectation value of a product of Wilson loops 
corresponding to a general link consisting of loops from this fibration. We will discuss this 
more below. 



2.2 Matter 

Next we would like to add matter to the theory. The matter will come in chiral multiplets, 
each of which consists of a complex scalar 0, a fermion ip, which is a 2-component complex 
spinor, and an auxilliary complex scalar F. 

The gauge-coupled action for a chiral multiplet in a representation R of the gauge 
group G is described in in the case of flat Minkowski space. It is straightforward to 
modify this for S^, giving: 

Sm = J d^xy/g (^Di,ct>^D^'4> + ^(t>^(t> + i^^p'ip + F'^ F - (t>^a'^(f) + cp^Dcj) - tp^faip + i</.UV - «V'"^A(/> 

(2.14) 

where, for example, Tp^atp is a shorthand for: 



V'tV^lT,)^^, (2.15) 

where a, b are indices in R, a is an index of the Lie algebra, and (T^)^ are the generators 
of in the representation R. Also, L*^ is a derivative that is covariant with respect to both 
the gauge group and the metric on S^, and we assume the various color indices have been 
contracted in a gauge invariant way. Note that the second term in ( |2.14| ), which arises 
from the conformal coupling of scalars to the curvature of S"^, gives the matter scalars a 



mass. A similar mass term will appear in the localizing term in ( 3.36 ). 



This action is classically invariant under the following superconformal symmetries: 

S(j) = rj^ip 

= {-i^D^cf) - ia4>) e - f 7^(V;,e)</' + 
5^^^ = r/t {i-ffD^(j)^ + i(T</.t) + i(/.t(v^?7t)^M + e^^t 

6F = {-i-i^D^il) + iXcj) + ia^)) 
SF^ = {iDf,^p''J^' - iX^(p^ + iaip^) r]* 



(2.16) 



Here e and i] must be a Killing spinors, satisfying ( |2.4| ). 

In order to perform the localization, it turns out that the theory must be super- 
conformal on the quantum level. This is because the supersymmetry 6 that we will use 
for localization, together with its conjugate and the Lorentz group, generate the entire 
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superconformal algebra. Thus any hermitian action invariant under 5 is necessarily super- 
conformal. 

This fact determines which superpotentials are allowed. In the absence of a superpo- 
tential, the combined Chern-Simons-matter system is superconformal on the quantum level 
1 13]. This follows from the nonrenormalization of the Chern-Simons couplings (except for 
finite shifts) together with the standard nonrenormalization theorem for the F-terms. An 
arbitrary quartic superpotential preserves superconformal invariance on the classical level, 
but in the quantum theory superconformal invariance is destroyed in general. Indeed, if 
the fields have anomalous dimensions, the scaling dimension of the quartic superpotential 
is not equal to 2, and the superpotential perturbation is not marginal. However, for spe- 
cial values of the superpotential couplings it may happen that the theory has enhanced 
supersymmetry which requires the anomalous dimensions to vanish [13|. This is the case 
for = 4 theories of Gaiotto and Witten, the Af = 6 ABJM theory, and the Af = 8 BLG 
theory. We will see later that the path-integral localizes to configurations where all matter 
fields vanish, so the precise choice of the superpotential will not matter, provided it ensures 
superconformal invariance on the quantum level. 

One particular example from the class of superconformal Chern-Simons matter theories 
is the ABJM theory Q. Here the gauge group is U{N) x U{N), with the Chern-Simons 
action for the two factors appearing at levels k and —k. The matter comes in two copies of 
the bifundamental representation (A'^, N), and two more in (iV, A^). There is also a quartic 
superpotential which ensures the supersymmetry is enhanced from J\f = 2 to M = 6. 



3. Localization 
3.1 Gauge Sector 

In this section we will closely follow |Q , in which supersymmetric Wilson loops were studied 
in = 2 and J\f = 4 super Yang-Mills theory in 4 dimensions, and their expectation values 
were computed by essentially the same localization method we use here. We will start by 
considering pure Chern-Simons theory, with no matter. We will discuss how the addition 
of matter affects the computation in section 3.4. 

The idea of the localization is as follows. We start by picking a single supersymmetry 
6 which preserves the operator we are interested in. We then deform the action by adding 
a term: 

t6V = MTr' (((^A)tA) (3.1) 

Here "Tr'" is some positive definite inner product on the Lie algebra.^ We assume this term 
is itself supersymmetric, which amounts to saying that V is invariant under the bosonic 
symmetry 5^. Then the standard argument shows that the addition of this term to the 
action does not affect the expectation value of any (5-invariant observable. 

^We distinguish it from the trace in the original action, which is not necessarily positive definite (eg, in 
ABJM, where it has a different sign for the two gauge groups). 
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We pick the term in (|3.1|) to deform the action because its bosonic part, t{6X)^6X, is 
positive definite. Thus we can take t to be very large, and the dominant contribution to 
the path integral will come from the region of field space where this term vanishes, which is 
precisely where 6X = 0. In the limit of large t the theory becomes free, so we can compute 
things easily, knowing the results we get are independent of t and thus apply at t = 0, 
where they represent the quantities we are interested in. 

Returning to the case at hand, let us fix a supersymmetric Wilson loop W along some 
great circle on S"^. This will be the operator we want to localize. 

We start by defining the supersymmetry we will be working with. Let e be the unique 
left-invariant spinor which preserves the Wilson loop, normalized so that e^e = 1, and let 
r] = 0. By "5" we will mean the infinitesimal supersymmetry variation corresponding to 
this choice of parameters. 

For the gauge multiplet, the transformations can be read off^ from ( |2.2| ) and ( p. 3D using 
r/ = and V^e = f 7^e. We find: 

6A^ = -§At7^e 
5a = -iAte 

6D =-i(D^At)7^e + iAte + i[At,a]e (3.2) 
5X = {-l^'^-'F^,- D + i-f>'D^a-a)e 
5X^ =0 

It is clear that 5^ = on the bosonic fields, and therefore on the fermions as well. 
Thus the (5-exact term ( |3.lD is trivially supersymmetric. As shown in appendix A, for the 
supersymmetry in question it evaluates to: 



SV = Tr' Qf^'^^F^. + D^aD'^a + {D + af + iX'^Y'^^.X + i[X\a]X - ^AU^ (3.3) 

Note that, unlike the matter scalars, there is no mass term for the scalar a. The mass 
term for a will arise not from the conformal coupling to the curvature of 5'^ but from the 
supersymmetric Chern-Simons action. 

Next we need to determine where the theory localizes to. The vanishing of 5X requires: 



= ( - -D + i^D^a -a]e (3.4) 



This implies the following two conditions: 

1 



(3.5) 

D = -a 
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It is straightforward to show that, on S^, the only solution to the first equation is to take 
F^u = and a = ao =constant, and then the second equation implies D = —Uo- 

Thus the theory localizes to the space of constant a, with D = —a and all other fields 
vanishing. One can also see this from the explicit form of 6V above. 

In the limit of large t, the exact result for the path integral becomes equal to the saddle 
point approximation: 

Z = j daoe'-^'^''^Zl,^Jao] (3.6) 

Here we integrate the contributions e^'''^'^"^ from the saddle points, which are labeled by 
cjo) together with the determinant factor ■2^i_;oop [co] coming from quadratic fluctuations of 
the fields about each saddle point. For the partition function, the classical contribution is: 

= -iJd^x^2Ti{ao^) (3.7) 
= - WTr (ao^) 

where we have used the fact that the volume of S*^ is 27r^. The Wilson loop gives an 
additional factor of: 



trijfe/'^"'^) =trjj(e2'^'^°) (3.8) 



Thus we have: 



Z = I duoew {-AiTT^Tv (^7,2)) Zf_i^Jao] 



(3.9) 



This is essentially the same result as was found in Q for M = 2 and = 4 SYM 
theory in four dimensions. For A/" = 4, it was shown that Zi^ioop = 1, so that the resulting 
matrix model is Gaussian, while for = 2 it was something more complicated, involving 
the Barnes G-function. In the next section we will find that, in our case, Zi_ioop is not 1, 
but it is still something relatively tractable. 

Before moving on, we should mention that, to be precise, we should really be localizing 
the gauge-fixed theory. That is, we should have started by introducing ghost fields c, c, 
and a Lagrange multiplier b. We would then have the standard BRST transformations 6b, 
and, continuing to follow [Q], define a new fermionic symmetry: 

5' = 6 + 6b (3.10) 
One can check that 5''^ = 0. We also would modify V to: 

V Tr' ({6X^)X) + cVAf, (3.11) 
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We would then localize with respect to 6' rather than 5. The 6' variation of the new 
V has four contributions: from 5 and 5b each hitting the two terms in V . For the first 
term, only 5 would contribute, since 5b is a gauge transformation and the term is gauge 
invariant, so the total contribution would just be the 5-exact term we found above. For 
the second term, the 5b variation would give us the usual gauge-fixing term: 

cV;,D^c + 6VM^ (3.12) 
Then we still need to worry about the remaining term: 

5{cV^'A^) (3.13) 

But if we define 5c = 0, then we are only left with some term multiplying c, which can be 
absorbed into the definition of c. 

In other words, we have shown that we can proceed by starting with the action in ( p.3| ) 
and gauge fixing it the usual way, and then computing the 1-loop determinant. We turn 
to that calculation now 

3.2 1-Loop Determinant 

In this section we compute the 1-loop determinant coming from quadratic fluctuations of 
the fields about the saddle points we found in the last section. After introducing ghosts, 
(|3.3|) becomes: 



S = t j ^d^xTi' (^^F'"'F^^ + D^aD''a + {D + af + ^X^]pX + i[X\a]X-^X^X + ^f,cD^'c + bV''Af,^ 

(3.14) 

We will be interested in the large t limit, so we rescale the fields to eliminate the t out 
front: 



D^-ao+^^D' (3.15) 



Here ^ represents all fields other than a and D, and we have treated these fields differently 
because they have zero modes. Also, a' represents the non-zero mode part of a, and 
similarly for D' . Taking t to be large then allows us to keep only the quadratic terms in 
the action: 



S= J Vad^xTv' ^^Fa^^'Fa^u + [A^. ao? + d^ad'^a + (D' + 

iX^yX + i[X\ao]X - ^AU + d^Wc + bVA^^ 



(3.16) 
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where Fa^u = dfj,Ai, — d^Afj_. The integral over D' can be performed immediately, and 
eliminates the squared term. The b integral is also easy, giving a delta function constraint 
which imposes Lorentz gauge. We find: 



(3.17) 

This is a free theory, and we would like to compute its 1-loop determinant. A similar 
calculation is done in [16|, and we can proceed similarly. First we separate the gauge field 
into a divergenceless and pure divergence part: 



A^ = d^<)) + B^ (3.18) 

where V^B^ = 0. Then the delta function constraint becomes 5(— V^(/)), and so we can 
integrate over (j) using the delta function, picking up a jacobian factor of det(— V^)^^/^. 
The integral over a' gives the same factor, while the integral over the ghosts contributes 
a factor of det(— V^). These all cancel (and in any case, are Uo-independent), and we are 
left with: 



S = j ^d^xTr' (^-B''V^B^ + [B^,aof + iX^fX + i[X\ao]X-^X^X^ (3.19) 

Now if we go back to ( |3.9D for a moment, we see that since the action is gauge invariant, 
the integrand is invariant under the adjoint action of the group. Thus we can replace the 
integral over the entire Lie algebra with an integral over some chosen Cartan subalgebra. 
This introduces a Vandermonde determinant in the measure, leaving, eg, for the partition 
function: 

Z = Jdal^Yl a{af^ exp {-AiTr^Tv {a')) Zl.^Ja] (3.20) 

where a runs over the roots of G, and a runs over the Cartan subalgebra. Thus we only 
need to know Zf_^^^p[ao] for cTq in the Cartan, and so from now on we will assume ao is in 
the Cartan. Now let us decompose as: 

B^ = ^ Bfj_°'Xa + 

a 

where Xq, are representatives of the root spaces of G, normalized so Tr' (XaXi^) = Sa/s, 
and a runs over the roots of G. Here is the component of along the Cartan, but 
this part of B^ will only contribute a Co-independent factor to the 1— loop determinant, so 
we will ignore it. Then we can write: 

[ao,B,,] = Y,a{^o)B^"X^ 

a 
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We can do something similar for A. Plugging this into the action, we can now write it 
in terms of ordinary (as opposed to matrix valued) vectors and spinors: 



S = j Vad^x [sf^a (-V2 + a{aof) B^^ + At, {i^ + ia{ao) - ^ \^ (3.21) 

From we know that the eigenvalues of the vector laplacian acting on divergenceless 
vector fields are (£ + 1)^, where £ = 1,2, and they occur with degeneracy 2l{l + 2). Thus 
the bosonic part of the determinant is: 

oo 

det(bosons) = H IK^^ + + "(^o)^)^^^^^^^ (3.22) 

a 1=1 

For the gaugino, we note that on 5^, eigenvalues of iy are ±(£ + ^) with degeneracy 
l{i + 1), where I runs over the positive integers.^ Thus the fermion determinant is: 

ia{ao)) (3.23) 

And so the total 1-loop determinant is: 

^l-loopi^o] =[[[[ 



a e=o ((£ + l)2+a(ao)^)^(^+2) 



a il (£ + ia(cTo))(^-i)(^+i)(£ + 1 - ia(ao))^(^+2) 



(3.24) 



We see there is partial cancellation between the numerator and the denominator, and 
this becomes: 

Because the eigenvalues of a matrix in the adjoint representation come in positive- 
negative pairs, only the even part in sigma contributes. We can isolate this by looking 
at: 



CO 

UU^^' + ^i^of? (3.26) 



a 1=1 



a(cJo)^ 2 



\l=l / a e=i 



This can be seen easily using the results of section 3.4, specifically as a special case of (3.53) 
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Ignoring the infinite constant, this product can be done exactly, and is simply: 



To summarize, we have shown: 



-f-..M-ni^«^j (3.-) 

where a runs over the roots of G. Plugging this into (|3.20| ), we see the denominator cancels 
against the Vandermonde determinant. Introducing the notation: 

detR/(a) = n/(p(«)) (3.29) 

p 

where R is some representation, and the product runs over its weights p (which, in the case 
of the adjoint representation, are just the roots of the algebra), we are left with: 



Z = J daexp (— 4i7r^Tr (a^)) det^if^ sinh(7ra) 

<:W > = ^ } ^ I daexp (-AiTr'^Ti (a^)) tiR (e^''") det^d sinhfyra) 
z dim it J \ / / \ / 



(3.30) 



The first line reproduces the result of [p, ^ that the Chern-Simons partition function 
can be obtained from a matrix model. Here it has been derived using the supersymmetry 
that the Chern-Simons theory possesses. 

3.3 U{N) Chern Simons Theory 

For a concrete example, we will look at the case where G = U{N). Then we can take the 
Cartan as the set of diagonal matrices, setting a = diag(Ai, Aat). The roots of G are 
labeled by integers i / j, and have: 

aij{a) = Xi- Xj (3.31) 

Also, as mentioned earlier, we take Tr as ^ times the trace in the fundamental rep- 
resentation, and we also take the Wilson loop in the fundamental representation. Then 



( |3.30D becomes (up to a sign): 



Z 



I (lle-'''-^'"dxA nsinh2(^(A, - A,)) 

\ i / i<j 



<3 

(3.32) 



TZ j (ll'""'^''''^') (e'''*'+-+e''*")n™''^('(^'-*j" 

\ i / i<j 
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To interpret this result, note that without matter, we can integrate out the auxiUary 



fields trivially. The integral over D in (2.1) imposes cr = 0, and so we see from (2.5) that, 
in the case of pure Chern-Simons theory, the super symmetric Wilson loop we have been 
considering is just the ordinary Wilson loop operator. Thus the second line of ( p. 32 ) gives 



a simple way of computing the Wilson loop expectation value in this theory. 

Both of the integrals above are just sums of Gaussian integrals, and it is straightfor- 
ward, though tedious, to evaluate them exactly, as shown in appendix B. The result for 
the Wilson loop expectation value is: 

sin (^) 

<W>=e-— . ^ l> (3.33) 

reproducing the known result up to an overall phase. 

This phase comes from the framing of the loop, which can be seen as follows. As 
mentioned earlier, the supersymmetry we are using preserves a family of Wilson loops 
forming a Hopf fibration of the sphere. The framing of a Wilson loop is essentially the 
choice of a nearby loop so that point-splitting regularization may be performed, and for 
this procedure to be compatible with supersymmetry, this loop must come from the Hopf 
fibration, and therefore have linking number —1 with the Wilson loop. 

One simple extension of this calculation is to a link of Wilson loops. If the loops in 
such a link come from the Hopf fibration, then they preserve the same supersymmetry 
5. Then it is easy to see that the expectation value of this operator is given simply by 
inserting more factors of tr/je^'^" in the matrix model, or equivalently, taking the trace in 
the product representation. This property of the Chern-Simons invariant of the Hopf link 
can also be shown by topological means. 

3.4 Matter Sector 

Next we carry out the localization procedure in the matter sector. Rather than treat the 
case of multiple chiral multiplets in various representations, we will consider a single chiral 
multiplet in some possibly reducible representation R. The action of the supersymmetry 5 
on the matter fields is as follows: 



(3.34) 



with all other variations vanishing. Here we have assumed the rescaling ( ^.15 ) has been 



done on the gauge multiplet. Since we will be taking t to be very large in a moment, this 
means we can ignore all coupling to the gauge sector, except that through Uq- 



^Thanks to Lev Rozansky for discussions on this point. 



- 14 - 



To localize the matter sector, we add a term similar to the one we used in the gauge 
sector: 



5V^ = t6 i^{5i;)H + V't(5^t)tj (3.35) 
In Appendix A it is shown that this equals: 

SVm = d^(t>^d''(t> + i(t>^v''df,^ + (t>^ao^(t> + ^(t>^(t> + F^F + ij^ (^J7 - tao + (^^-^^"j V (3-36) 

where v'^ = e^j^e. As before, this term is positive definite, and vanishes on the following 
field configurations: 

0=5-0 = {-i^i^dfj_(l) - iooct) + \(t))e 

(3.37) 

= 5-0"^ = e^F^ 

The second equation implies -F = 0. With a little work, one can check that the first 
implies (p must also be zero. Rather than show this directly, we will see in a moment, when 
we evaluate the 1-loop determinant, that the operator acting on (p in ( p.36 ) has no zero 
modes, which leads to the same conclusion. 

Thus there is no classical contribution coming from the matter sector, and its only 
influence is through its effect on the one- loop determinant. Since there is no interaction 
between the matter and gauge flelds at quadratic order, except that throught Uo, the 
determinant f actor izes, and we can write: 

Z = j daexp (— 4i7r^Tr a^) det^(i(sinh(7ra))Z]^^QQp[a] (3.38) 
We will compute this extra factor in the next section. 
3.5 1-Loop Determinant - Matter Sector 

For the scalar field, we see from ( |3.36| ) that the operator we need to diagonalize is: 



Dbos = + iv^'d^ + - + ^o^ (3.39) 



1 
4 

Here "do" is actually a matrix representing ao in the representation R. As we did for the 
gauge multiplet, we can decompose this representation into its weight spaces. Namely, 
if Bp is a representative of the weight space corresponding to the weight p, satisfying 
< Cp^Cpi >= 5ppi (where < .,. > is some gauge- invariant way of contracting the relevant 
color indices), then we write: 



(3.40) 



The total 1-loop determinant will be the product of the one coming from each term in this 
sum, which are all acting on ordinary (not matrix-valued) scalars. 
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It will be most convenient to use a pair of orthonormal frames, one left-invariant and 
one right-invariant under the action of SU (2) (thinking of 5^ as SU (2) and letting it act 
on itself). We will call these T and r*. Then we can take v = l^. It is straightforward to 
show that the laplacian can be expressed in terms of these fields as: 

V' = = E(^*)' (3.41) 

i i 

where we think of the vector fields as differential operators on the space of scalar fields. 
Thus the each term in (|3.40| ) can be written as: 

Dbos = -kl' + i/^ + ^ + p{aof (3.42) 
Also, using the fact that the vectors satisfy the algebra: 

%lj] = -2eijklk (3.43) 

we see that if we define new operators Lj = — these satisfy the SU{2) algebra. In terms 
of these operators, the operator acting on the scalars becomes: 

= 4L,L^ - + 1 + p{ao)^ (3.44) 



and so computing its eigenvalues reduces to a familiar problem from quantum mechanics. 

e 

"2 



On a spin- 1 representation, the determinant of the operator can be written as: 



det,/2(Aos)= n (^(^ + 2) - 2m + ^ + p(ao)2^ (3.45) 

m=-| 

It can be shown that the scalar fields on 5^ decompose into the irreps (|, |) under 
the action of the left- and right-acting SU {2)'s, so the total determinant will be a product 
of the above expression over all non-negative integers i, each raised to the power of the 
degeneracy, which is ^ + 1 owing to the right-acting SU{2). 

Next we consider the fermions. After decomposing them into weights as with the 
bosons, the operator we need to diagonalize is: 

Dferm = +\i ^P^S'o) (3.46) 

If we use the li as our vielbein, the covariant derivative acting on spinors can be written 

as: 



Then the Dirac operator is: 



= + -7m (3.47) 



iV = - \ (3-48) 
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We should be careful to distinguish 7*/j, which is a differential operator, from j) = 7^, which 
is just a matrix. Thus the operator acting on the fermions becomes: 

Dferm = ifk - 1 + " ^P(^o) (3.49) 

Or, if we define = 7V2, which satisfy the SU{2) algebra, and plug in the Lf 

= -AS'Li + S^-1- ip{ao) (3.50) 

So the problem reduces to computing spin-orbit coupling. Unfortunately, we cannot 
proceed the standard way since does not commute with the total angular momentum 
J = L + S, so we are forced to compute the determinant manually. Let: 

= 2aL -3 + 2^83 + J (3.51) 

Note that this operator commutes with both and J^, so its eigenvectors all have 
the form: 

V = a\-,m > \-t > +b\-,m + l > \i> (3.52) 
Letting O act on these vectors, it is straightforward to compute: 



dete/20 = {a^ + P + j){a^-^ + j) J] (-^(^ + !)«' - (2m + l)a/3 - 07 - + 7') 

m=-| 

(3.53) 

Plugging in the relevant values, a = —2,P = ^,7 = — 1 — ip{ao), we get: 

det^/2(D/e.„,) = l + ip{<^oW+l + ip{<7o)) n (^(^ + 2) - 2m + ^ + p{aof^ 

m=-| 

(3.54) 

We note this is almost equal to the scalar determinant, except for the extra terms out 
front and the missing the m = | factor in the product. Taking this into account, we can 
write: 

deti/2Dferm ^ , ^ ^ (^ + ^ + ip(ao)) (^ + | + ip(ao)) 
dete/2Dbos e{£ + 2)-£+l + p{aof 

But the denominator factors as: 



e{e + 2)-£ + -+ piaoY = {£+ -f + piaoY = + 2 + ipi^o)){£ ip{(To)) (3.56) 

So we see one of these factors cancels one of factors in the numerator, and we are left with: 
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Then the fuh determinant is given by taking the product of this over ah i, keeping 
track of the I + 1 degeneracy at each level. Letting n = £ + 1, and also taking the product 
over eigenvalues p{cro), we get: 

detDferm TT TT + I + ^P('^o ^ 



deti;^, p^=i V" - 2 - / 
To deal with this quantity, we start by looking at its log: 



log {Zr^ioopi'^o]) = X] £ ^ f'' + ^ + ipi^^o)] - log ("n - ^ - ^ (3.59) 

p n=l ^ ^ ^ ^ ^ ^ 

To isolate the (To-dependent piece, we take a derivative with respect to cJo- Letting cij be 
the components of do in some basis of the Cartan, and defining pi = p{ei), we have: 



_d_ 

da. 



p n=l - (2 

(3.60) 

+ i/9(o-o))' 



.2 



„2 _ /'i 

'^n=l / p n=l 



1+^/5(0-0))^ 



The second term converges, and we will compute it in a moment. But the first term 
seems to present a problem, since it contributes a o"o-dependent divergence. Specifically, if 
we integrate back in the o"o dependence, we find: 



lo: 



g (^r-zoopko]) = ( E 1 J 2. E p^^^o) + ... 

\n=l J p 

(3.61) 

= (^E^^ 2iixR{ao) + ... 

If the gauge group G contains no U{1) factors, or if the representation is self-conjugate, 
then tri{((To) will vanish and so this term will not contribute. However, even if it does not 
vanish, it is not hard to see that we can absorb this term into a renormalization of the 
Fayet-Iliopoulos term in the action. 

Thus we only need to worry about the second term in ( p.60|) . The sum over n can be 
done in closed form, and we find: 



da. 



log (zr_,„„p[ao]) = i ^ (^1 + vri Q + ip{aS^ tanh(7rp(^7o))^ (3.62) 
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In general, this cannot be integrated in terms of elementary functions to get back 
log ^-Z^j^ioop^ • The problem is the imaginary part. If we assume that i? is a self-conjugate 
representation (ie, as a set, {p} = {—p}), then we see that the imaginary part cancels in 
the above sum, and we are left with: 

p 

1 d (3-63) 
= "2 0^^^°^ (cosh(7r/j(ao))) 

So that, up to a (To-independent constant: 

ZT-ioop[< = n (cosh(7rp(cTo)))-i/2 ^3_g4) 
p 

We will focus on the case of self-conjugate representations here, but we note in passing 
that we can still apply the localization to more general representations, at the cost of an 
extra phase in the 1— loop determinant: 

p 

where /(x) satisfies ^ = 1 ~ xtanhx, and can be expressed through the polylogarithm 
function. 

Thus we have shown, in the notation of (3.29): 



^™ ioopko] = detij (cosh(7rcro)) (3.66) 

where this result only applies when the matter appears in a self-conjugate representation. 
A simple example of such a representation is one of the form R = S (B S* , in which case 



( 3.66 ) becomes: 



^T-ioopWo] = dets (sech(7rcro)) (3.67) 

We will restrict to such representation for the remainder of this paper. 

To summarize, we have shown that the partition function of the super symmetric Chern- 
Simons theory with chiral multiplets in representations Ri, R^, R2, R2--- localizes to the 
following matrix integral: 

Z = y daexp (— 4i7r^Tr a^) det^^ (sinh(7ra)) det^^ (sech(7ro)) det^j (sech(7ra)) ... (3.68) 
while the expectation value of the supersymmetric Wilson loop operator (p.5|) is: 



<W>= 

dim 



— — — : / daexp (— 4z7r^Tr a^) tr/j(e^™)detAd (sinh(7ra)) det/j^ (sech(7ra)) det/jj (sech(7ra)) 
iRZ J 

(3.69) 
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4. Results for ABJM Theory 



As an explicit example of the results of the last section, we will look at ABJM theory. This 
will provide a check of our calculation, as the Wilson loop we are interested in has already 
been studied perturbatively elsewhere . 

To be precise, the Wilson loop considered in those papers was the following: 

W = ^Tr Pexp {J dr (i^^x^ + X^Xslil)) ) (4.1) 

Here X"^ are the scalar fields of the theory, of which there are four, and Xa are their 
adjoints. is a constant hermitian matrix which can be taken as diag(l, 1, —1, —1). In 
those papers, it was shown that this choice of M renders the Wilson loop 1/6 BPS, ie, it 
preserves one real supersymmetry and one superconformal symmetry (they were working 
in flat Minkowski space). While this is the same as for the Wilson loop we have been 
considering, it is not obvious they are the same operator. 

However, as shown in |17] (specifically, equation (4.11)), the quantity appearing in 



the second term of (4.1) is precisely what we get for a after integrating it out. Thus this 
is the same as the operator we have been considering, only written in a form where the 
supersymmetry is less manifest. 

As mentioned earlier, the matter content of ABJM theory is two chiral multiplets in 
the bifundamantal (A^, N) representation, and two more in the dual (iV, N) representation. 
Writing a = diag(Ai, Aat, Ai, Xn), the roots run over all z, j from 1 to A, and satisfy: 

Pij (a) =Xi-Xj 

(4.2) 

For the ABJM theory the scalar product which we denoted TV is 

A(tr-tr), (4.3) 

where tr and tr are traces in the fundamental representation of the two U{N) factors of 
the gauge group, and k is an integer. Plugging this into (|3.68| ) we find that the partition 
function localizes to the following matrix integral: 



f . \ ni<, (smh2(^(A, - A,)) sinh2(vr(A; - A,))) 

Z= [ T[e'^^^^^ -^^ ^dXidXi ^ (4.4) 

J V¥ 7 n„cosh2(vr(A.-A,)) 

The Wilson loop is given by inserting ^ • e^'^'^^ as before. This is no longer a simple 
Gaussian integral, and we were not able to obtain an exact result for general N . 
In order to perform a perturbative calculation, we can rewrite this as: 
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xexp I E.<, I 21og ( > ) +21og I ^ I I -2E„log(cosh2 

(4.5) 

where we have defined A(A) = HiKjl'^j ~ ^i)^ ^"^^ ^ — ~^ — 2vrii, where t = ^ is the 
't Hooft coupUng. Also we have ignored overall constants. The second line is regular as 
A — )■ 0, and so can be expanded as a power series in A. The first line is a product of two 
Gaussian measures: 

'"n(e"^^''fiAi) A(A)2 (4.6) 

i 

The expectation value of powers of A in this model can be computed exactly using 
orthogonal polynomials as shown in Clearly odd powers of A have vanishing expec- 
tation value, and it is not hard to see that < >= 0{a^). Thus if we are interested in 
computing expectation values in a small a expansion, it is sufficient to consider only the 
first few terms in the small A expansion of the second line above. 

Carrying out this procedure, we find: 

1 1 ^- I \ 2 1 /, , 4 \ 3 



<»r>=i + -„--(^i + ^j„^--(^i + ^j„. + ... (4.7) 



From what we saw in ( 3.33| ), we would expect this to be the result for a framing of 



— 1. To compare with the perturbative field theory calculation, where trivial framing is 
assumed, we will need to determine and remove this phase. But since the linear term is 
due entirely to this phase, we can accomplish this simply by multiplying above result by 
This gives the result for trivial framing: 

Plugging in a = 2Trit, we arrive at the final answer: 

/5 1 \ TT^N^ (\ 1 \ ivr^iVS 



In the large N limit, the second order term agrees with the result of IC, 11|. 
5. Discussion 



The main result of this paper, ( 3.68 ), can be applied to any of the theories discussed in the 



introduction. It may be possible to compute some of the resulting matrix integrals exactly, 
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as we did for pure Chern-Simons theory. If not, this result would still provide a much more 
efficient way of performing perturbative calculations. 

Another possibility is that these matrix integrals simplify in the large A^-limit. Since 
this is the limit we are usually interested in for the AdS/CFT correspondence, exact results 
(as a function of the 't Hooft coupling) in this limit would be almost as good as exact results 
for finite N. One approach one might take for these large N calculations is the saddle point 
method described in |12|. 



It should also be possible to extend this method to calculate the partition function 
and supersymmetric Wilson loops on more general manifolds. In particular, the partition 
function for Chern-Simons theory on Seifert manifolds was also shown to reduce to a matrix 
model in [P, ^, and one would expect this could be shown using a localization calculation 
similar to the one we performed on S^. 

A. Appendix: Q-exact Terms 

In this appendix we derive the following results. The gauge sector localization term is: 

dV = Tr' (^^F'^'F^, + D^aD^'a + {D + af + iAt7'^V^A + i[X\a]X - ^A+A^ (A.l) 
while the matter sector term is: 

'1 + ' 



(A.2) 



5Vra = {d^<P^d^'<j) + iv>'<P^d^(l) + ^^ao^(l) + ^(/^V + F^F + V'^ (^J7 - iOo + 
A.l Gauge Sector Calculation 

We have (ignoring the trace for notational convenience): 

W = <5(((^A)U) (A.3) 

From (3.2), the supersymmetry transformation b has the following action on the gauge 
multiplet fields: 

bA^ = -iAt7^e 
5a = -iAte 

5D =-i(D^At)7A^e + |Ate + i[At,a]e (A.4) 
5\ ={-\Y''F^^-D + ^-^^'D^a-a)e 
(5At =0 

The bosonic part of the Q-exact term is: 
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5Vbos = (5A)t(,5A) 

/I \ / 1 \ ^^-^^ 

= (-jl'^'F^^- -D- i^D^a - a) {-^r^F^, -D + ij^'D^a - aj e 

Note that the signs are such that the F + Da terms do not mix with the D + a terms, 
and we are left with: 

SVbos = \f''''F^, + D^aDi'a + {D + af (A.6) 



The fermionic part is: 
We compute: 

= d^6A^ - d„A^ + i[A^, 5A^] - i[A^, 5A^] 



Df,dA^ - DJAf, 



-i(L>^(At7,e)-Z?,(At7^e)) 

-i (p^At)7,e + |(At7,7^e) - {D^X^h^e - ^iX^^^jue)) 
-|(I)^At)7,e + |(D.At)7^e + i(At7,^e) 



and: 



6{Di,a) =S{di,a + i[A^,a]) 
= D^6a + i[5A^,a] 
= L»^(-^Ate) + i[-iAt7^e,a] 
= -l(D^At)e- iAt7^e+ i[At,a]7;.e 

Plugging this in gives: 



(A.7) 



(A.8) 



(A.9) 
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(^~{D,X^)l^e+\\h+'-[\\a]^ + (A.IO) 

These term naturally fall into 3 groups: Those with covariant derivatives, of which 
there are 4, those with cr's, of which there are 2, and those with neither, of which there are 
4. The first group gives (combing the first two using the antisymmetry of 7'^'^): 

,t (_!_^,^^^D,\^)^,e + '-{D,X^)re + \r{D,\^)e'^ A (A.ll) 
Now we will use the following Fierz identity (for anticommuting spinors): 

(r?i^'?2)(??3^?74) = ~{r]i^m){'n2}m) - ^{mhiJ.m){v3h'^V2) (A.12) 

This gives: 

= -k{- i(eh'^'^7.e)(i^^At)A - i(et^M-^^^^e)(i^^At)7PA- 

+i(ete)(Z)^At)7MA + i{e^^,e){D^\^)^^j-^ X- (A.13) 

meU){Df,X^h^X + i{eh.e){D^X^)Yj''x) 

Using 'j'^'^'ji, = 27^ and 7'^'^7p7i/ = — 2(5p, we see the first, fourth, and sixth terms 
cancel while the rest combine to give simply: 

-zV^A^7'*A = iAt7'^V^A (A.14) 
Next we look at the second group in SVferm: those involving a. These are: 

e^(^-'-[X\a]e-'-[X\a]j,e^X (A.15) 
But now we can apply the Fierz identity the other way to get: 

= ^(e^e)[A^(7]A = ^[A^a]A (A.16) 
Finally, the last group gives: 

et Q^M^(At^^^,) _ 1 (Ate) - l7''(At7^e) + l(Ate)) A (A.17) 
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Using 'j^'^ = ie^'^P^p, we see the first term combines with tlie third, and we get: 



(i7nAt7^6) + l(Ate)) A 
-\ (i(eh^7Me)AU + 1 (et^A«^^^^e)At7-A + i(ete)AtA + \{e^ ^,,)\^ Y (A.18) 



-iAtA 



Putting this all together, the total 5-exact piece becomes: 

5V = ^F^^'F^^ + D^aD^'a + {D + af + iAt7^V^A + i[\\a]\ - ^AU (A.19) 
A. 2 Matter Sector Calculation 

In the matter sector, the (5— exact term splits into two pieces: 

5Vm = 5Vml + SVm2 



For the first piece, the transformations we will need are (from ( |3.34| ) ): 



(A.20) 



(A.21) 



6iP = (^-iYd^4) - iao(j) + |0)e 

where we have assumed we have already localized the gauge fields, so that the only coupling 
to the gauge multiplet is through ctq. The bosonic part is: 



SVmibos = {S06^p 

= et(i7'^9^(/.t + i(^^ao + ^(t>^){-ij^'df,4> - iao^ + i0)e 

+|(et7Me)((5^0t)<^_^t(a^^)) 

(A.22) 

where = e^'^^e. Using the constancy of cTo and the divergencelessness of v^, the fifth 
term is a total derivative, and the second and last terms can be integrated by parts to give: 

= a^</.ta^</) - ie^'^p{VpVp)<i)'^d^(t) + (t)'^ao^<j) + - it;^</,ta^0 (A.23) 

Finally, using V pVy = e^ypV^ , we find: 



^Vmuos = 5^0^ + iv^<t)^d^(t> + cp^o^cp + i(/.t0 (A.24) 
Next we compute the fermionic part: 

= et(i7'^9^(V^te) + ia^ii^^e) + i(V^e))^ (A.25) 
Fierz rearranging, using ( |A.12| ): 



(A.26) 

Integrating by parts: 

(A.27) 

But using: 

(V^7;,)7'^7" = ^^.upvPYl^ = e^.p^''ie'^^'^7<x = -^iv^'^^ (A.28) 
We see the second term combines with the fifth term, and we get: 



So that the final answer can be written as: 

SVml ferm = V'^ (^)^ + 1 " ^^o) V' (A.30) 

Note this is degenerate. However, we have not yet locahzed F. This will come from 
the second term: 



5Vm2 = S(i>^{6i>^y) (A.31) 
The relevant supersymmetries here are: 
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6F =€^i-i-f''V^,^ + iao'ip) 

(A.32) 

5V;t = e^Ft 



The bosonic part is simply: 



(A.33) 



while the fermionic part is: 

SVm2ferm = i^^'SFe* 

(A.34) 

Noting following relations: 



e^e* = ete = 1 



(A.35) 



The above becomes: 

6Vm2 ferm = (^^) " (^-36) 

SO that it combines with the fermionic part of the other term to give: 

SVmferm = V'^ (^f - iUo + ) V' (A.37) 

Putting this all together, we have: 



38) 



B. Appendix B: Explicit Computation of Partition function and Wilson 
loop for U{N) Chern Simons Theory 

In this section we compute the following integrals: 

Z = J (He-^'-^^'dxA nsinh2(7r(A, - A,)) (B.l) 

\ i / i<j 
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\ i I i<j 

We will use the following identity: 

Yl 2sinh(xi - Xj) = e-(^+i)(^i+-+^^) J^i-l^e^^J "^^^^ (B.3) 

l<i<j<N a 

where a runs over all permutations of N elements. This is essentially a special case of the 
Weyl character formula. Using this identity, we get for the partition function: 



f (Yle-^'-^'dxA I]sinh2(7r(A, - A,)) 

\ i ) i<j 



<3 

(B.4) 

2 



2N{N 

Putting a = ikir, P = 27r, this becomes: 

0'1,IT2 \ i / 

This is just a sum of ordinary Gaussian integrals (provided we add an appropriate 
small imaginary part to k to ensure convergence), and performing the integrals gives: 

0-1 ,0-2 \ i / 

(B.6) 

= \ /vr\iV/2 ^ gEi(^i«+^2«-(iV+l))2 

2iV(JV-i) \a) l^^ ^> 

cri,<T2 

Consider the sum in the exponent: 



N 



Y,Mi) + a2{i)-{N + l)f 



i=l 



(B.7) 



TV 



{ai{if + a2{if + {N + if - 2{N + l)ai{i) - 2{N + l)a2{i) + 2ai{i)a2{i)) 



i=l 

All but the last term are just sums of powers of the integers from 1 to A'', and it becomes: 



= -^N{N + l)(iV + 2) + 2^ ai(f)a2(i) (B.8) 

i 

This leaves: 

0"l,cr2 

For fixed ai, the sum over a2 can be performed using ( [B.3D , yielding: 



27V(7V-i) Va 



1 ^vr^A^/2^_^^(^^^^(^^2)^(_i).i^|^(7V+i)K(i)+...+<^i(iV)) -Q 2sinh(^(ai(i)-ai(j))) 



(B.IO) 

Finally, since the product in odd is a, the sum simply introduces a factor of A^!, and we 
get: 



i<«<i<A^ 

(B.ll) 

m=l ^ ^ 



2Af(7V 

Plugging in for a and /3, we arrive at (ignoring an overall /^-independent phase): 

- - n (=) (B-) 

m=l 

Next we turn to the Wilson loop. Since the integrand for the partition function is 
symmetric under permutations of the Aj, this will be given by the following integral: 

r (Yle'^^^'^dxA e^^^'llsmh\7r{Xi - A,)) (B.13) 

\ i J i<j 

Starting at ( |B.5| ), we see all that changes is that, in the exponent, there is an extra 
term linear in Aj with coefficient /3(5ii. Thus the sum in the exponent is modified to: 



^ 1 

^(fTi(i)+a2(i)+5ii-(A^+l))' = --iV(iV+l)(A^+2)+2^ai(i)a2(z)+2(cTi(l)+<J2(l)-(iV+l))+l 

1=1 I 

(B.14) 

Plugging this in gives: 



-29- 



2N{N-i) U; Z^l e e 

(B.15) 

We have written it this way so that we can still perform the sum over a2 using ( |B.3D , which 
gives: 



2A'(A^-l)/2 



a 



(B.16) 



xj](-ire^'^« n sinh(^(a(i) + 5H-^7(i)) 

Note we do not need to include 8ij in the sinh since j > 1. Now, unless cj(1) = N, there 
will be some j such that the argument of the sinh is zero. Thus we only need to consider 
those terms with cr(l) = A^. Then this becomes: 



^ (_l)-e^A^/'jJsinh('^(iV+l-a(.)))) f n smh(^(cT«-a(i))) 
o-eSjv-i Vi=2 ^ ^ / \2<i<j<N ^ ^ 



\m=l ^ ^ / o-eSjv-i \2<i<j<N ^ ^ 



m=2 ^ ^ / \ l<i<j<7V-l 



sinh ( 1^ 



^ sinh ( ^ 



(B.17) 



Thus we see the normalized Wilson loop expectation value is: 



„2 sinh (^^^ 

<W>= e^^ ^^""^^^ 

sinh 1^ 



Or, plugging in the correct values of a and /3: 



-30- 



„■ (nN\ 

-TTiN Sin —f— 

< W >=e^ 

sin (f ) 



(B.19) 
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